Quasi-ID atomic gases across wide and narrow confinement-induced-resonances 
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We study quasi-one-dimensional atomic gases across wide and narrow confinement-induced- 
resonances (CIR). We show from Virial expansion that by tuning the magnetic field, the repul- 
sive scattering branch initially prepared at low fields can continuously go across CIR without decay; 
instead, the decay occurs when approaching the non-interacting limit. The interaction properties es- 
sentially rely on the resonance width of CIR. Universal thermodynamics holds for scattering branch 
right at wide CIR, but is smeared out in narrow CIR due to strong energy-dependence of coupling 
strength. In wide and narrow CIR, the interaction energy of scattering branch shows different types 
of strong asymmetry when approaching the decay from opposite sides of magnetic field. Finally we 
discuss the stability of repulsive branch for a repulsively interacting Fermi gas in different trapped 
geometries at low temperatures. 
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I. INTRODUCTION 

Quasi-onc-dimcnsional(lD) atomic gases across scat- 
tering resonances can be realized in laboratories by uti- 
lizing the confinement-induced-resonance(CIR)[l[. By 
initially preparing the system at high or low magnetic 
field and sweeping the field properly, the quasi-lD sys- 
tem can evolve on the attractive branch with molecules [21] 
or on the repulsive scattering branch that is free of 
molccules[3, ^. Quasi-ID atomic gases have many fasci- 
nating properties that are very different from those of 3D 
gases. For instance, in quasi-lD a two-body bound state 
exists for an arbitrary s-wave scattering length a^ [l|, I4I . 
The three-body recombination rate of ID bosons is effi- 
ciently suppressed in the Tonks-Girardeau(TG) regime 
with strong repulsion^, Q- A long-lived metastable 
quantum phase in the super-TG regime [y, 0] has been 
realized in the scattering branch of a ID bosonic system 
with strong attraction [J]. 

Apart from dimensionality, the resonance width is an- 
other important ingredient affecting many-body proper- 
ties. Take (3D) Feshbach resonance(FR) for example. 
In wide FR, where the width is much larger than the 
typical energy scale £*[8], the system exhibits univer- 
sal thermodynamics (UT) right at resonance where as 
diverges |9l [lO|. UT means that the thermodynamic po- 
tential is a universal function of the temperature and den- 
sity, regardless of any detail of inter-particle interactions. 
However, in narrow FR, where the width is much smaller 
than £* , the universality is not evident due to consider- 
able effective-range effect [ll|. Another interesting prop- 
erty in narrow FR is that the interaction energy shows 
strong asymmetry when approaching resonance from dif- 
ferent sides [131 ■ as recently been observed in ^Li Fermi 
gas[lj|. Considering the facts that the quasi-lD geome- 
try is reduced from 3D by confinement and the ID res- 
onance originates from the 3D s-wave interaction, it is 
natural to expect the effective-range effect in 3D will also 
influence the interaction properties of quasi-lD system. 

In this work, from the analysis of two-body solu- 
tions and high-temperature Virial expansions, we study 



the scattering property and thermodynamics of quasi-lD 
atomic gases across CIR. We shall show that the quasi-lD 
geometry greatly modifies the stability of repulsive scat- 
tering branch compared to 3D case. Similar to FR, CIR 
can also be classified as wide or narrow according to the 
resonance width. We find drastically different thermody- 
namic properties between wide and narrow CIR due to 
the energy dependence of coupling strength. The stabil- 
ity of repulsive branch in other trapped geometries, such 
as isotropic 3D trap or anisotropic quasi-low-dimensional 
trapped system, is also discussed in combination with re- 
cent developments on cold Fermi gases in the laborato- 
ries. We use h = kB = 1 throughout the paper. 

The paper is organized as follows. In Sec. II, we intro- 
duce effective quasi-lD scattering, from which the wide 
and narrow CIR are defined and the corresponding ther- 
modynamics are presented. In Sec. Ill we carry out 
high-T Virial expansion for effectively ID system, and 
present a detailed study on the stability and thermody- 
namic properties of the repulsive branch. An extensive 
discussion of the stability of repulsive branch for a cold 
Fermi gas in various trapped geometries is given in Sec. 
IV. Finally we summarize the paper in Sec. V. 



II. EFFECTIVE ONE-DIMENSIONAL 
SCATTERING AND THERMODYNAMICS 

The Schrodinger equation for the relative motion of 
two atoms moving in quasi-lD is 



FoV'(r) + ^^^^^<5(r)|-[rV'(r) 
m or 



,0 = E^Pir), (1) 



here the non-interacting part is Hq — — V^/ni + 
muj'^{x^ + 2;^)/4, r — {x,y,z) and r = |r|; In the 
pseudopotential part, we use the energy- dependent s- 
wave scatteriiig length obtained from a renormalization 
procedure [l4l [l5j. 



as{E) = abg{l + 



W 



E/S^L -{B- Bo) 



(2) 



as (E) physically describe both wide and narrow FR, with 
background scattering length ai,g, magnetic field B, reso- 
nance position Bq, width W, and magnetic moment dif- 
ference Sfi between the atom and closed molecular state. 
The reduced quasi-lD scattering from 3D s-wave in- 
teraction has been solved by Olshanii et al[l|. For low- 
energy scattering with E = uj±_ + k'^/m and k'^/m <C 
2a; _L, the wavefunction at large inter-particle distance 
is frozen at the lowest transverse mode, and its even- 
parity part is phase-shifted as is ^cvcn(r) ~ exp[— (x^ -I- 
y^)/{2a\)] cos(fc|z| -I- Sk), with a± — ^/2/{mujjJ and 



cotSk 



ka\ 



'as{E) 



Co + o(- 



-)], (3) 



where Co = 1.4603. Hereafter we neglect the small 
correction from last term in Eqj31 5k in turn deter- 
mines an ID energy-dependent coupling strength, g{E) = 
2fctan5fe/m with E = k'^/m, as 



g{E)=gbg{l 



WiD 



E/S^i -(B- Bid) 



), 



(4) 



where gbg = 2-/uj±abg, Wm = ^W, Bid = Bq - 
1)W — uj±/6fi, with 7 = (1 — Coabg/a±)~^{see also |15l|). 
EqU explicitly shows all realistic parameters describing 
CIR, namely, the background coupling g^g, resonance po- 
sition Bid and width Wid- Near CIR {B ^ Bid) and 
for E ^ SfiWiD, one can construct an effective-range 
model to formulate ID interaction. 



9iE) 



1 

9lD 



= TT^ roE, 



(5) 



with giD the zero-energy coupling strength and tq = 
— l/{mabgSfiW) the effective range characterizing E- 
dependence in as(£') [lg|. To this end, Eas. (|4l5|) show the 
reduced effective-range effect (or E-dependence of cou- 
pling strength) from 3D to quasi-lD system. 

In the tight transverse confinements and low atomic 
density(7i) limit, na± <?C 1, we consider an effective ID 
system with interaction given by Eq|4l Generally, the 
pressure takes the form 

P (0 \^^f^ , Sii{B-Bid) 5^iW,d Ebg 
P = ^(2m/i)2 jr(-, { , , — =i}), 

n ^J. tJ. tJ- 

(6) 
where /i is the chemical potential, E^g = Tag?, and J- is 
a dimensionless function. 

For wide CIR, (5/iW^i£)(> 2uj±) > n^/m, the E- 
dependence in EqH] and EqlS] is negligible, and the in- 
teraction parameters in {...} of Eq|B]can be replaced by 
a single giD- The pressure is then reduced to 



P = f,{2mtJ.)iTi-,-^). 



(7) 



At wide CIR{giD = oo), P is just a function of T and u 
(or T and n) indicating UT for the scattering branch|Tr]. 
Particularly at T = 0, UT can be established by noting 



that the bosons and spin-1/2 fermions with infinite repul- 
sion are fully fermionalized, with the ener gy i dentical to 
that of an ideal single-species Fermi sea[3,[l^|. However, 
at narrow CIR(i? = Bid), EqlS] still essentially relies on 
other interaction parameters {Wid, gtg) and thus UT is 
absent. More explicitly, UT can be identified by Virial 
expansions at high temperatures. 



III. HIGH-TEMPERATURE VIRIAL 
EXPANSION 

At high temperatures, n? /ra <^ T <^ 2ix)j^, we c arry 
out virial expansions on the effectively ID system TS]. 
The pressure can be expanded using the small fugacity 
z = e^/^ as P = a^ E„>i ^n-^". where A = ^27r/(mr) 
is the thermal wavelength, a is 1 for spinless boson and 2 
for equal mixture of spin-1/2 fermions. Compared with 
non-interacting case (with superscript "0"), 



P = p(") 



T 



Y.^K-1^^)-^ 



(8) 



n>2 



Here the difference, hn — hn, characterizes the interaction 
effect to the n-body cluster, which is generally a function 
^^|MB_Biii)^ Sj^^ :^}. For wide CIR, &„-6i°^ only 
depends on one single parameter ll{XmgiD), which is 
free of parameter at giD = oo for any order of Virial ex- 
pansion and leads to UT according to EqlHl This also jus- 
tifies us in examining UT within the second-order Virial 
expansion. Consideration of higher-order expansions will 
not change the conclusion, except for a negligible correc- 
tion (of higher order in z or nX) to the thermodynamic 
quantities. 

Due to the interaction effect, the second Virial coef- 
ficient, A62 = (62 — ^2 )/v^, can be written as A62 = 

^Je~^'/"^ — e^^i l^\ (here I is the energy level for rel- 
ative motion of two atoms). Given P(T,[i) in EqlH it 
is straightforward to obtain the density n = dP/d^i and 
entropy density s — dP/dT , and finally energy densities, 
8 = fin -\- Ts — P, for spinless bosons(b) and spin-1/2 
fermions (f) as 



£f = 



n.T 

2 

nT 



nX , 
^57j(-l + 2e„,) 



o{{n\f) 



1 

l + -^(l + 2e„,,) + o((nA)^) 



%\ 



25/2' 



with dimensionless interaction energy 

.aA62 



tint = -A62 



2T- 



dT 



(9) 
(10) 

(11) 



In the following we derive A62 in strictly ID by enu- 
merating the energy levels of two interacting particles 
in a tube([— L/2, iy/2]). For simplicity, we first consider 
the scattering branch without inclusion of any bound 



state. The discretized wavevector(fc > 0) is determined 
by boundary condition 

kiL/2 + 5i = {I + \/2)ti (/ = 0,1,...). (12) 

By comparing to non-interacting fc^ where 5i — 0, we 
obtain 

A62 = ^[exp(-fcf/(mr)) - cM-k^i'^VimT))], (13) 



which can be transformed to an integral 
2/(mT) /g°° dfcfc(5fce-'='/('"^) and further to 



as 



A6? 



dke 



"^ dk 



(14) 



Note that to obtain EqUHwe extrapolate Si^q to dk=o in 
the thermodynamic limit, and set i5fc=o — — 7r/2 consider- 
ing di=o < as well as EqlS) When considering a bound 
state (occupy / = 0), the lowest available I for scattering 
state should he I = 1. This implies Sk=o is up-shifted by 
TT as revealed by Levinson's theorem. In this case, 



A6! 



bd 



^\E,\/T 



dke 



, dSk 
'dk' 



(15) 



To this end A62 is obtained for both repulsive scattering 
branch fEa lTil) and attractive branch (EqfT5)). Remark- 
ably, comparing with 3D 9], A62 in ID has an additional 
term (—1/2) resulted from zero-energy phase shift and 
the unique scattering property of ID system. Eqs. (|14|15p 
are consistent with results obtained from Bethe-ansatz 
solution[2(j| and analyses of real-space wave functions j2l|. 
For quasi-lD system, the Virial expansions are carried 
out by setting k — 2/a± as an upper limit of integrals 
in Eqs. ()14ll5p . In the rest of this section, we shall mainly 
focus on the scattering branch(cf Egfl^ which might ex- 
hibit UT as discussed above. 
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FIG. 1: (Color online). Upper panel: Schematic plot of zero- 
energy coupling strength gm across wide (Al) and narrow 
(A2) CIR. The labels (a-e) correspond to B < Bio(a);_B -^ 
Bid -0+{b);B ^ Bm +0+(c);B > BiD{d);B > Bmie). 
Lower panel: Phase shift Sk versus k across wide (Bl) and 
narrow (B2) ID resonances, with each label (a,b,c,d,e) corre- 
sponding to specific giD as marked in (Al) and (A2). 



here x = 2^/2n/{m\giD\)^)', sgn() is sign function and 
erf() is error function. In the weak coupling limit {x — >■ 
00), we obtain 



A6r = -l/iV^x), e- = 2/(V^x) 



(17) 



for scattering branch at gm — > 0+ (corresponding to 
solid lines in small B-field in Fig. 2); and 



A^^'^ = l/(V^a;), 



hd 



= -2/(V^x) 



(18) 



for attractive branch at gm -^ 0~[2^ (dashed lines in 
large B-field in Fig. 2). In the strong coupling limit [x — ?> 
0), we obtain 



With (5/iWi_D(> 2ijj±) > n^/m, we replace the E- 
depcndcnt g{E) by a constant gm- giD is schematically 
plotted in Fig.l(Al), giving the phase shift [SfS) of scat- 
tering branch in Fig.l(Bl). Here we have excluded the 
existence of bound state for any B-field, thus 5k all start 
from — 7r/2 at /c = 0. By increasing B across CIR(from 
" a" to " e" ) , the amplitude of 5k at finite k gradually be- 
comes enhanced, implying more repulsive energies in the 
system. In particular, 5k is uniformly — 7r/2 for all k right 
at CIR(between "b" and "c"), leading to universal values 
of A6|'^ and ef^^ as shown below. 

For strictly ID system with constant giD, Eqs. p^fTS)) 
can be analytically solved, for example. 



A6? 



— 2 — '^^p(;;2)[i-^^/(-)]' (16) 



A5? 



1 1 , x-^ 

— ± — ^(x 

2 20F^ 2 

1 x^ 

2^27^' 



(19) 
(20) 



the universal values at x = 0, — Afel'^ = ef^^ = 1/2, are 
direct consequences of fc— independent phase shift (— 7r/2) 
as mentioned above. 

In Fig. 2, we plot A62 and tint for two-species ^Li 
fermions across wide CIR. For scattering branch, we see 
that all curves of A6|'^(or ef^J at different T intersect 
at a single point in A&2'^(or ef^J-B plane [23|, demon- 
strating the UT of scattering branch right at wide CIR. 
The scattering system at strongly repulsive side of CIR 
can smoothly evolve to strongly attractive side with even 
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2: (Color online). A62(a) and ei„t(b) for two- 
species Li fermions across wide CIR at T{iiK) = 
1.5(dark black), 3(niedium purple), 6(liglit pink). We con- 
sider FR at Bo = 834. IG with width W = -300G[2il. 
Transverse confinements are generated by optical lattices 
with lattice space ul = SOOriTTi and depth Vo ~ 25Er 
{Er = ^(-^f), giving uji_ = (27r)300KHZ. CIR occurs at 
Bid = -Bo -^152. 2G with Wm = -147.9G, which satisfies 
S^iWiD S> 2uj± 2> T. Solid and dashed lines are respectively 
for scattering(Ea fT4l) and attractive (Eg [TS]) branch. Dash-dot 
lines denote universal values — Abl'^ = ef^t — 1/2 at CIR. 
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higher energy. This is consistent with previous theoret- 
ical predictions of sTG phasep, 0] and its recent exper- 
imental realization in bosonic gas[4|. Virial expansion 
also shows that the scattering branch will achieve the 
strongest repulsion as g\u -^ 0~ at large B-field, with 
— A&l'^, e|^( — T' 1. All above properties can be clearly seen 
by tracing any individual energy level of two scattering 
atoms in a tube, as shown in Fig. 3(a). 

Here we remark on the stability of scattering branch. 
In the framework of two-body clusters in Virial expan- 
sion, the decay of scattering branch manifest itself in the 
discontinuity of thermodynamic quantities, due to the re- 
labeling of scattering states when the underlying bound 
state converts to the lowest scattering state. This is why 
in 3D the decay occurs right at FR where the bound 
state converts to scattering state at a^ = oo^9|. In ID, 
however, the conversion is at g\D — instead of at res- 
onance, and therefore the scattering branch can extend 
far away from CIR until approaching zero coupling limit. 
A more comprehensive discussion of the stability of scat- 
tering branch in other trapped geometries will be given 
in Section IV. 

At the end of this subsection we briefly discuss the 
second-order virial expansion in a ID harmonic trap, 
which can be carried out given the two-body spectrum 
under coupling strength Eq|4l In particular, at wide CIR 
with g\Y) = -|-oo, the spectrum is Ei = (21 + 5/2)ujz com- 
pared with EJ"'' ^ {21 + 1/2)lUz; this gives Ab^2%ap = 
— 1/(2a/2) compared with —1/2 in homogenous case. In 
fact, based on local density approximation(as used in 3D 



FIG. 3: (Color online). Two-body energy levels in the 
center-of-mass frame for quasi- ID system confined in a 
tube([— L/2, L/2]) across wide (a) and narrow (b) CIR. The 
orange and black dashed lines denote 7r/2 and phase shift 
(corresponding to g{E) = oo and 0). Eq = {2it/L)^ /m. 
The dotted lines denotes non-interacting energy levels with 

i5(«'/-Eo = (^ + i)^ / = o,i,2.... 



trapped system in Ref . |24| ) . we have 
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(21) 



for the scattering branch right at wide CIR. This shows 
a more rapid convergence of virial expansions in trapped 
ID system than in homogeneous case. 



B. Narrow CIR 

With S^iWiD ^ m? /m{<^ wj,), we take the full form 
of g{E) (EqU due to the strong E-dependence. As- 
sume a positive background ai,g, we give the schematic 
plot of giD in Fig.l(A2) and 6k for scattering branch in 
Fig.l(B2). In Fig. 4, we show A62 and ei„t for ^^Rb sys- 
tem across extremely narrow CIR. 

Compared with wide CIR case, the scattering branch 
in narrow CIR shows many distinct properties. 

First, 5k is no longer universal at CIR; instead, it sen- 
sitively depends on k and is quite small at finite-fc (see 
"b" and "c" in Fig.l(B2)). This can be_attributed to the 
strong E-dependence in EqU that g{E) is far off reso- 
nance at finite E even its zero-energy value gio — ^ 00. 



Accordingly, as plotted in Fig. 3(b) the two-body levels 
at finite energies are just shifted by a small amount al- 
though the lowest level is shifted half-way. As a result, 
there is no UT at narrow CIR, and the scattering branch 
is generally with very weak repulsion even close to CIR 
(see also Fig. 4). 

Second, shortly beyond CIR, the scattering branch 
goes through a decay a.t B — Bid + Wm as manifested 
by discontinuous Afe^'' and ef^^ there(see fig. 4). This is 
exactly the place where gio evolves from 0~ to 0"*" and 
the bound state transforms to scattering state. 

Third, after the decay, i.e., B > Bid + Wm and 
giD > 0, we see from Fig. 1(B2) that Sk will complete 
a continuous change from — 7r/2 to nearly tt within an 
energy window AE k, B ~ Bid- Due to the large tt shift 
for all energies larger than Ai?, we see an large and neg- 
ative ef^t in Fig. 4(b) despite of positive giD- Similarly 
to narrow FR in 3D 12], we expect the negative e,f,'^( in 
extremely narrow CIR will extend to much larger B-field, 
until 5ii{B — Bid) approaches the typical energy scale of 
the systemfS|. 

On the whole, the scattering branch in narrow CIR has 
weak repulsion(ef^j — > 0+) or strong attraction(ef^( -> 
(— l)"*") when B approaches the decay position {B = 
Bid + Wxd) from the small or large field side. The 
asymmetry here differs from that in wide CIR, where 
ef^( approaches 1~ or 0+ respectively. It is also helpful 
to compare these features in quasi- ID with those in 3D 
system [3, Il2|- For wide FR in 3D, the amplitudes of in- 
teraction energies are symmetric for the system evolving 
in different branches and approaching FR from different 
sides (9|, which is in contrast with what we find in quasi- 
ID system across wide CIR. For narrow FR, the inter- 
action effect are greatly suppressed for repulsive branch 
but greatly enhanced for attractive branch[i2|, the same 
features as revealed above in quasi- ID system across a 
narrow CIR. In cold atoms experiments, all these features 
in quasi-lD system can be detected using the technique 
of rf spectroscopy, as has been successfully applied to a 
quasi-2D Fermi gasf20l and a 3D Fermi gas across narrow 
FRpj]. 



IV. STABILITY OF REPULSIVE FERMI GASES 
IN TRAPPED GEOMETRIES 

Recently, the metastable repulsive branch of atomic 
gases has attracted lots of research interests, in the con- 
text of the experiment by MIT group on itinerant fer- 
romagnetism for a repulsively interacting Fermi gas[27J. 
The same group later claimed the absence of itinerant 
ferromagnetism from the measurement of spin suscepti- 
bility, and attributed this to the instability of repulsive 
branch against the molecule formation for a 3D Fermi 
gas near Feshbach resonance [2g|. Similarly, the instabil- 
ity of a repulsive Fermi gas has also been observed in the 
quasi-2D Fermi gasQ, but at negative Og side. 

There have also been quite a few theoretical studies 
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FIG. 4: (Color online). A62(a) and ei„t(b) for 

**^Rb bosons across narrow CIR at T{fj,K) = 
l(dark black), 2(medium purple), 3(light pink). We con- 

sider FR at Bo = 406. 2G with width W = OAmG^. 
The optical lattice is same as that for ®Li in Fig.2, giving 
cjx = (27r)80KHZ. CIR occurs at Bid = So - 27.1mG 
with WiD = 0.5mG, which satisfies Sfj.WiD -^ T <^ 2luj_. 
Decreasing B across Bid + Wid , the scattering branch(solid 
lines) continuously evolves to attractive branch(dashed lines) 
with a bound state emerging at threshold. 



as to why the repulsive Fermi gas in a 3D homogenous 
system is unstable close to Feshbach resonance |29l - [3l| . 
For instance, the instability has been attributed to the 
shifted resonance in the background of a Fermi sea|29|. 
the pairin g in stability dominating over ferromagnetism 
instability '3(j|, or the vanishing zero- momentum molecule 
due to Pauli-blocking effect of Fermi-sea atoms ^3 UJ- All 
these studies can lead to the same conclusion at low tem- 
peratures, i.e., the 3D Fermi gas becomes unstable at 
the place where Og is comparable to the inter-particle 
distance {1/kp), or equivalently, the two-body binding 
energy {Eb ^ l/mal) is comparable to the Fermi energy 
{Ep ^ kp/m). The scattering branch can only be stable 
when Ef) > Ep, where the deep molecule can not be ab- 
sorbed by the Fermi sea atoms ^29|-l31|. In other words, in 
this parameter regime the existence of a deep two-body 
bound state effectively stabilizes a many-body system at 
the metastable repulsive branch. Since the physics be- 
hind this criterion does not depend on any detail of the 
dimension or trapped geometry, it should be equally ap- 
plicable to other cases besides the homogeneous 3D sys- 
tem. In the following, we will use this criterion to study 
the stability of repulsive Fermi gas in various trapped 
geometries at low temperatures. 

Typically we consider three different types of trapping 
potentials, namely, the isotropic or nearly isotropic 3D 
trap (lJx '^^ Wj, ~ Wz), the extremely anisotropic quasi-2D 
(wz ^ iLjx,ujy) and quasi- ID trap (w^ '^ uix,ujy). To fa- 
cilitate the discussion, we consider the system across wide 
resonance (with single interaction parameter), while the 
extension to narrow resonance should be straightforward. 

In a trapped system, a two-body bound state is always 



sumjorted no matter how weak the attractive interaction 
is[E[32,|33]. For a 3D isotropic trap (uj^ ~ Wy ~ w^ ~ w), 
however, it should be noted that the two-body binding 
energy Ei, at a^ < side is less than the order of to, 
i.e., the level spacing of all scattering states[32]. As a 
result, in the thermodynamic limit with atom number 
N :^ 1, Ei,{< w) at as < side is negligible com- 
pared with the Fermi energy Ep ~ Ncu. The system 
is then expected to behave similarly to the homogeneous 
case, in a sense that the repulsive branch is stable only 
with positive Ug where the bound state is visibly deep 
{Eb ~ Nuj). This is consistent with what have been ob- 
served in MIT experiments [27|, l28|. On the contrary, for 
anisotropic quasi- 2D or quasi- ID trap, the energy spac- 
ing of scattering state is generally of the order of trapping 
frequency of the shallow confinement, while the binding 
energy can be of order of trappin g fr equency of the tight 
confinement even at a., < side[l|,|33|. For example, for a 
quasi-2D trapped system at a^ = oo, Ei, > ujz '^ uj,j;,ujy, 
and the existence of deep molecule would be possible 
to stabilize the repulsive branch at a^ = oo as long as 
Eb > Ep ~ Nu}x,y In this case, the stable scattering 
branch can even extend to negative a^ side, as shown in 
the experiment with a quasi-2D Fermi gas|26|. The same 
conclusion can be drawn in the quasi-lD trapped case 
[ujz <S^ ujxi'^y), which is also consistent with the high- 
temperature result presented in the last section. 

In short summary of this section, at low temperatures, 
the stability of repulsive Fermi gas in a trapped geometry 
relies not only on the existence of two-body bound state, 
but more importantly, on the value of its binding energy 
compared with the typical energy scale of a many-body 
system. In other word, here the two-body physics should 
be evaluated in a many-body background. Generally, 
the repulsive branch in quasi-low-dimensional systems is 
expected to be more stable than that in an isotropic 3D 
system. Therefore the low-dimensional system provides 
us a more favorable platform to realize possible itinerant 
ferromagnetism in repulsively interacting Fermi gases. 



V. SUMMARY 

In this paper, we have studied the quasi-lD atomic 
2:ases across wide and narrow CIR. Our main results are 



summarized as follows. 

First, from high-temperature Virial expansions we ob- 
tain the following: 

(1) By tuning the magnetic field across CIR, the re- 
pulsive scattering branch of quasi- ID system can evolve 
continuously across CIR, from gio = +oo to gip = — oo 
side. 

(2) Universal thermodynamics are identified for the re- 
pulsive scattering branch right at wide CIR, but is found 
to be washed away at narrow CIR by the strong energy- 
dependence of coupling strength. 

(3) The decay of quasi-lD repulsive branch occurs 
when gi£) — > 0. The interaction energy shows differ- 
ent types of strong asymmetry between wide and narrow 
CIR, when approaching the decay position from opposite 
sides of magnetic field. 

Moreover, the second-order virial expansion presented 
in this paper also serves as a benchmark for testing future 
experiment on ID atomic gases. 

Second, we have discussed the stability of repulsive 
branch for a repulsively interacting Fermi gas at low tem- 
peratures in different trapped geometries. By evaluating 
the two-body bound state in the presence of a Fermi sea, 
we conclude that the system can generally be more stable 
in the quasi- low-dimensional trapped system than in a 3D 
isotropic trap. This should shed light on the current ex- 
periments seeking for ferromagnetism in more stable and 
strongly interacting Fermi gases in low dimensions. 
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